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ABSTRACT 

This  paper  describes  the  evaluation  of  the  NRL  Limited 
Area  Dynamical  Weather  Prediction  Model  in  simulating  coastal 
atmospheric  mesoscale  phenomena  along  the  west  coast  of  the 
United  States  during  the  period  0000  UTC  02  May  -  1200  UTC  03 
May  1990.  A  graphical  comparison  technique  was  used.  Model 
output  was  compared  horizontally  with  large-scale  analyses, 
station  data,  cross-section  analyses,  and  vertical  profiles  at 
specific  locations.  The  model  successfully  simulated  the  wind 
and  temperature  fields,  but  failed  to  accurately  replicate 
moisture  and  height  fields. 


in 


TABLE  OF  CONTENTS 


I.  INTRODUCTION   1 

II.  MODEL  DESCRIPTION  AND  DATA  SET 4 

A.  BASIC  MODEL  DESCRIPTION  4 

B.  INPUT  DATA 9 

III.  METEOROLOGICAL  SITUATION  (02-03  MAY  1990)   ....  13 

A.  SYNOPTIC  SCALE  ANALYSIS  13 

B.  MESOSCALE  ANALYSES   18 

IV.  MODEL  PERFORMANCE  FOR  SYNOPTIC  SCALE  PHENOMENA   .  .  24 

A.  GENERAL 24 

B.  ANALYSIS  AT  0000  UTC  02  MAY  1990   (MODEL  INPUT 
FIELDS) 26 

1.  Constant  Pressure  Surfaces  26 

a.  Surface 26 

b.  850  Mb  Constant  Pressure 28 

c.  700  Mb  Constant  Pressure 30 

d.  500  Mb  Constant  Pressure 31 

2.  CROSS -SECTION  ANALYSES  31 

C.  MODEL  OUTPUT  AT  1200Z  02  MAY  1990 34 


IV 


^wiyut  i  i\ivua  LlbKARY 
NAVAL  POSTGRADUATE  SCHOOI 
MONTEREY  CA  93943  5101 

1.  Constant  Pressure  Surfaces  34 

a.  Surface 34 

b.  850  Mb  Constant  Pressure 39 

c.  700  Mb  Constant  Pressure 41 

d.  500  Mb  Constant  Pressure 45 

2.  Cross- section  Analyses  45 

a.  MFR  to  NSI 45 

b.  NSI  to  ELY 50 

3.  Vertical  Profiles   55 

D.  MODEL  OUTPUT  AT  0000  UTC  03  MAY  1990 59 

E.  MODEL  OUTPUT  AT  1200  UTC  03  MAY  1990 60 

V.  MODEL  PERFORMANCE  FOR  MESOSCALE  PHENOMENA  62 

A.  GENERAL 62 

B.  LAND/SEA  BREEZES  AT  FORT  ORD,  MONTEREY,  CALIFOR- 
NIA    63 

C.  THE  SOUTHERLY  SURGE  ALONG  THE  CALIFORNIA  COAST  65 

D.  TOPOGRAPHICALLY  INDUCED  GRAVITY  WAVES  68 

E.  THE  CATALINA  EDDY 69 

VI.  CONCLUSIONS 75 

LIST  OF  REFERENCES 76 

INITIAL  DISTRIBUTION  LIST  79 


ACKNOWLEDGEMENTS 

I  would  like  to  extend  my  sincere  thanks  to  Professor 
Teddy  Holt  for  his  guidance  and  support  throughout  the  course 
of  this  study.  His  patience  was  truly  remarkable.  A  special 
thanks  goes  to  Professor  Wendell  Nuss  for  critical  review  of 
this  thesis.  Thanks  also  to  Jim  Cowie  and  Craig  Motell  for 
their  programming  assistance,  and  to  Russ  Schwanz  and  Donna 
Byrch  of  the  Naval  Postgraduate  School's  IDEA  Lab  who  often 
came  in  during  their  off -hours  and  kept  the  IDEA  Lab  running; 
their  support  was  instrumental  in  completing  this  work  on 
time.  Thanks  to  Dick  Lind  who  provided  the  NPS  profiler  data. 
Special  thanks  to  Chi -San  Liou  of  NRL -Monterey  for  providing 
the  initial  model  fields,  and  to  Simon  Chang  of  NRL  for  model 
development.  Most  important,  however,  thanks  to  Carolyn  -  my 
best  friend  and  lifelong  inspiration. 

This  research  was  funded  by  the  NPS  Direct  Funded  Merit 
Program.   Technical  sponsor  was  the  NRL-Monterey. 


VI 


I .   INTRODUCTION 

It  has  long  been  recognized  that  a  large  portion  of 
significant  weather  occurs  on  the  mesoscale  (Orlanski  1975) . 
The  recent  advent  of  large,  fast  computers  has  stimulated  the 
development  of  limited  area  or  regional  numerical  prediction 
models  capable  of  capturing  these  mesoscale  events.  One  such 
model,  the  Naval  Research  Laboratory  (NRL)  Limited  Area 
Dynamical  Weather  Prediction  Model,  has  developed  continuously 
since  its  initial  inception  in  the  1980 's  and  is  now  used 
confidently  to  diagnose  mesoscale  phenomena  and  their  interac- 
tion with  larger  scale  systems  (Holt  et  al .  1990;  Schulz 
1992) .  Although  the  NRL  regional  model  has  been  run  success- 
fully for  various  geographic  locations  and  domains  throughout 
the  world,  it  had  yet  to  be  used  for  the  western  United  States 
area. 

Within  the  western  U.S.  region,  coastal  and  topographic 
interaction  greatly  influence  the  structure  of  the  planetary 
boundary  layer.  Specific  west  coast  topographic/coastal 
mesoscale  studies  include  the  Catalina  Eddy  event  (Bosart 
1983;  Mass  and  Albright  1989)  with  its  associated  southerly 
surge  (Dorman  19  87)  and  windward  damming  of  winds  or  frontal 
features  in  the  Pacific  Northwest  (Mass  and  Ferber  1990) . 

Of  all  these  studies,  however,  few  prognostic  modeling 
studies  have  been  done.   Ulrickson  and  Mass  (1990)  performed 


numerical  simulations  of  a  Catalina  Eddy  using  the  three 
dimensional  Colorado  State  University  Mesoscale  Model  in  order 
to  investigate  synoptic  influences  on  pollutant  transport  in 
the  Los  Angeles  Basin.  Schultz  and  Warner  (19  82)  and 
Glendening  et  al .  (1986),  using  a  two  dimensional  mixed  layer 
model,  also  focused  on  the  transports  of  pollutants  in  that 
basin.  These  studies  were  confined  to  the  relatively  small 
domain  of  the  California  bight  and  did  not  look  beyond  for 
greater  spatial  distribution  of  topographic  and  coastal 
influenced  mesoscale  events. 

The  purpose  of  this  paper  is  to  validate  the  ability  of 
the  NRL  Limited  Area  Dynamical  Weather  Prediction  Model  to 
replicate  atmospheric  mesoscale  features  which  occurred  along 
the  west  coast  of  the  United  States  during  the  period  0000  UTC 
02  May  1990  to  1200  UTC  03  May  1990.  The  model  had  not  been 
used  previously  to  simulate  sea/land  breezes  and  mountain  or 
valley  winds  comparable  to  those  found  in  the  California 
region,  so  no  guidance  was  available  concerning  model  perfor- 
mance or  modifications  required.  Although  several  statistical 
techniques  have  been  developed  to  quantitatively  assess  model 
results  (Keyser  and  Anthes  1977;  Willmott  et  al .  1985) ,  a  less 
rigorous  technique  of  graphical  and  tabular  comparison  was 
used.  This  technique  was  chosen  for  several  reasons.  First, 
it  allowed  evaluation  of  predicted  fields  over  the  entire 
domain  rather  than  at  specific  observation  locations.  Second, 
analysis  fields   (which  might  have  served  as  observation 


fields)  were  based  on  a  synoptic  observation  network  which 
could  not  completely  resolve  mesoscale  phenomena.  Any 
conclusions  drawn  from  statistical  comparison  between  analysis 
and  forecast  fields  would  therefore  be  confined  to  the  larger 
scale. 

In  Chapter  2  of  this  paper  the  3-D  limited  area  model  is 
described  along  with  input  data  and  topography  fields  used  in 
the  experiment.  Chapter  3  describes  observed  synoptic  and 
mesoscale  features  which  occurred  during  the  period  02-03  May 
1990.  The  model's  performance  in  handling  large  and  mesoscale 
features  is  detailed  in  Chapter  4  and  5,  respectively. 
Conclusions  are  included  in  the  last  chapter. 


II.   MODEL  DESCRIPTION  AND  DATA  SET 

A.   BASIC  MODEL  DESCRIPTION 

The  Naval  Research  Laboratory  (NRL)  Limited  Area  Dynamical 
Weather  Prediction  Model  is  a  mesoscale,  baroclinic,  hydro- 
static three  dimensional  model.  Incorporating  processes 
affecting  both  small  and  larger  scale  phenomena,  it  is  appro- 
priate for  modeling  the  lower  troposphere  in  which  the  balance 
of  large-scale  motions  are  approximately  gradient.  A  brief 
description  of  the  model  is  provided  below;  more  extensive 
model  descriptions  are  provided  by  Madala  et  al .  (1987)  and 
Chang  et  al .  (19  89) . 

The  model's  seven  governing  primitive  equations  are  in 
surface  pressure  weighted  flux  form.  Because  of  it's  terrain 
following  characteristics,  sigma,  defined  as  the  ratio  of 
pressure  to  surface  pressure  (Phillips  1957) ,  is  used  as  the 
vertical  coordinate.  Twenty  three  sigma  levels  resolve  the 
vertical  as  indicated  in  Table  1.  Five  of  the  governing  equa- 
tions are  prognostic.  They  include  the  u-  and  v-  momentum 
equations,  the  thermodynamic  equation,  the  moisture  continuity 
equation,  and  the  surface  pressure  tendency  equation.  The 
remaining  two  equations,  the  hydrostatic  and  continuity  equa- 
tions, are  diagnostic.   These  equations  form  a  closed  system 


Table  I.   Model  Sigma  Levels. 


Model  Level 

Sigma 

Del  Sigma 

1 

0.05 

0.1 

2 

0.15 

0.1 

3 

0.25 

0.1 

4 

0.35 

0.1 

5 

0.45 

0.1 

6 

0.55 

0.1 

7 

0.64 

0.08 

8 

0.715 

0.07 

9 

0.78 

0.06 

10 

0.835 

0.05 

11 

0.88 

0.04 

12 

0.915 

0.03 

13 

0.94 

0.02 

14 

0.957 

0.014 

15 

0.969 

0.010 

16 

0.978 

0.008 

17 

0.985 

0.006 

18 

0.99 

0.004 

19 

0.9935 

0.003 

20 

0.996 

0.002 

21 

0.99775 

0.0015 

22 

0.999 

0.0010 

23 

0.99975 

0.0005 

for  the  seven  dependent  variables  u,  v,  T,  q,  p.,  0,  and  da/dt. 

The  equations  are  approximated  by  a  first -order  accurate 
finite  difference  scheme  in  space.  This  form  enhances 
geostrophic  adjustment  while  controlling  small-scale  gravity 
waves.  A  C-grid  (Arakawa  and  Lamb,  1977)  is  used  for  the  fi- 
nite differencing.  This  grid  best  simulates  geostrophic 
adjustment  while  conserving  integral  properties.  The  model 
has  103  X  91  horizontal  grid  points  with  uniform  1/6  degree 
resolution  in  latitude  and  longitude  over  the  domain  28°-43°N, 
130°-113°W.  Temperature,  geopotential ,  specific  humidity  and 
sigma  are  computed  at  mass  points  (i,j)  while  east -west 
velocity  is  computed  at  the  midpoint  of  mass  points  along  the 
x-axis  and  north- south  velocity  is  computed  at  midpoints  along 
the  y-axis.  Vertical  velocity  is  evaluated  at  half -levels  in 
the  vertical  between  sigma  surfaces. 

Time  integration  utilizes  a  split-explicit  scheme.  This 
method  allows  a  time  interval  equal  to  four  times  that  of  a 
conventional  leapfrog  scheme  (Chang  et  al .  1989)  by  applying 
varying  time  steps  for  different  modes  of  the  linearized 
governing  equations. 

It  is  generally  agreed  that  improper  treatment  at  the 
boundaries  of  a  regional  model  can  offset  any  advantage  gained 
by  higher  resolution.  A  temporal  relaxation  scheme  is  used  in 
which  the  values  within  five  grid  lengths  of  the  boundaries 
are  relaxed  toward  the  large-scale  analysis.    Boundary 


conditions  are  updated  every  12  hours  and  nudged  towards 
hourly  interpolations.  This  scheme  is  preferred  over  less 
capable  fixed,  time -dependent,  or  sponge  schemes  (Chang  et  al . 
1989  and  Holt  et  al .  1990). 

Parameterized  physics  in  the  model  include  convective 
and  nonconvective  precipitation,  dry  convective  adjustment  and 
a  parameterized  planetary  boundary  layer. 

A  modified  Kuo  scheme  (Kuo  1974)  is  used  to  parameterize 
convection  precipitation.  Convective  precipitation  occurs 
when  grid  scale  low- level  moisture  convergence  exists  in  a 
convectively  unstable  environment.  This  low- level  moisture 
convergence  either  increases  the  humidity  of  the  air  or  is 
condensed  and  precipitates  as  rain.  Per  Anthes  (1977) ,  the 
partitioning  of  latent  heating  and  moistening  is  determined  by 
the  column  mean  relative  humidity.  Vertical  distribution  of 
heating  is  determined  by  the  difference  of  temperature  between 
the  environment  and  the  pseudoadiobat . 

Nonconvective  precipitation  is  reached  when  saturation 
occurs  on  the  resolvable  scale.  The  Clausius-Clapeyron 
equation  is  used  to  compute  excess  moisture  and  isobaric 
heating.  Depending  on  the  level  at  which  saturation  occurs, 
excess  moisture  either  precipitates  into  the  lower  layers  and 
evaporates  or  falls  to  the  surface. 

In  order  to  neutralize  superadiobatic  lapse  rates,  dry 
convective  adjustment  parameterization  is  included  in  the 
model.   The  adjustment  can  take  place  over  several  layers 


above  the  planetary  boundary  layer  (PBL) .  This  adjustment  is 
activated  when  the  static  energy  of  a  layer  exceeds  that  of 
the  adjacent  higher  layer.  A  slightly  stable  lapse  rate 
results  while  total  static  energy  is  conserved. 

PBL  parameterization  uses  the  turbulent  kinetic  energy 
(TKE)  E  -  e  closure  scheme  (Gerber  et  al .  1989;  Holt  and  Raman 
19  88) .  The  surface  boundary  layer  is  parameterized  by  Monin- 
Obukov  similarity  (Businger  et  al .  1971).  The  large  gradi- 
ents in  topography  associated  with  the  Sierra  Nevada  precluded 
the  stable  calculation  of  surface  roughness  length  as  a 
function  of  terrain  heights;  therefore,  roughness  lengths 
were  assigned  a  constant  value  of  5  cm  everywhere  over  land. 
This  value  is  commonly  used  for  most  operational  numerical 
weather  prediction  models.  Over  water,  Charnock's  rela- 
tionship was  used.  In  addition,  the  model  incorporates  a  soil 
slab  model  (Blackadar,  19  76)  to  predict  ground  temperature 
based  upon  the  surface  energy  equation  (Chang  1979) .  Initial 
sea  surface  temperatures  (SST) ,  however,  were  held  constant  in 
time  throughout  the  model  integration. 

A  second- order  horizontal  diffusion  scheme  is  included  in 
the  model  to  suppress  computational  instability.  The  non- 
dimensional  diffusion  coefficient  (approximately  0.004)  is 
sufficiently  small  to  enable  the  model  to  remain  nearly 
undamped,  but  still  prevent  the  growth  of  nonlinear  insta- 
bility in  the  solution.  Horizontal  diffusion  of  temperature 
and  mixing  ratio  do  not  contain  diffusion  in  the  vertical, 
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rather,  vertical  diffusion  occurs  on  terrain  following  sigma 
surfaces  over  mountainous  topography. 

B.   INPUT  DATA 

The  basic  data  set  used  in  the  initialization  was  the 
archived  Fleet  Numerical  Oceanography  Center's  NOGAPS  (Navy 
Operational  Global  Atmospheric  Prediction  System)  2.5  degree 
global  analyses.  Initial  fields  were  sea  level  pressure 
(SLP) ,  sea  surface  temperature  (SST) ,  D- values,  u- components 
(U)  and  v- components  (V)  of  velocity,  temperature  (T)  ,  and 
vapor  pressure  (e)  .  These  input  fields  are  listed  in  Table  2. 
D-values  were  converted  to  geopotential  heights  using  the  NACA 
standard  atmosphere  relations  (Haltiner  and  Martin,  1972) . 

Prior  to  initialization,  NOGAPS  fields  were  compared  with 
the  National  Meteorological  Center  (NMC)  final  Global  Data 
Assimilation  System  (GDAS)  analysis.  The  only  difference 
noted  between  the  two  analyses  were  in  the  moisture  fields, 
particularly  at  lower  levels.  Although  satellite  imagery  was 
not  available,  upper  air  data  suggests  that  the  GDAS  moisture 
analyses  may  be  more  correct. 

In  order  to  initialize  the  model,  fields  were  interpolated 
both  horizontally  and  vertically.  Horizontal  interpolation 
was  by  cubic  polynomial  to  1/6  degree  resolution.  Bilinear 
interpolation  was  utilized  along  the  boundaries.  Fields  were 
then  vertically  interpolated  to  the  model's  23  sigma  levels. 
Temperature  was  interpolated  linearly  in  log  pressure.  Mixing 


Table  II.   NOGAPS  Gridded  Input  Fields 


Level  (MB) 

Parameters 

Surface 

SLP,  SST 

1000 

D- Value, 

U 

,  V 

e,  T 

925 

D- Value 

U 

,  V 

,  e,  T 

850 

D-Value, 

U 

V 

e,  T 

700 

D-Value 

U 

V 

e,  T 

500 

D- Value 

U 

,  V 

e,  T 

400 

D-Value 

U 

,    V 

e,  T 

300 

D-Value 

U 

,  V 

,  e,  T 

250 

D-Value 

U 

,    V 

T 

200 

D-Value 

U 

V 

T 

150 

D-Value 

U 

,  V 

T 

100 

D-Value 

U 

,    V 

,  T 

70 

D-Value 

U 

,  V 

r   T 

50 

D-Value 

U 

,    V 

,  T 

30 

D-Value 

U 

,  V 

r   T 

20 

D-Value 

U 

,  V 

,  T 

10 

D-Value 

,  u 

,  V 

r   T 

10 


ratio  was  interpolated  exponentially  in  pressure.   U  and  V 
velocity  components  were  interpolated  linearly  in  pressure. 

Topography  was  derived  from  the  U.S.  Navy  global  10 -minute 
elevation  data.  Five  point  horizontally- averaged  smoothing 
was  applied  to  the  elevation  data.  This  provided  a  near- 
realistic  terrain  field  for  the  model  validation  run.  Figures 
1  and  2  show  domain  wide  and  hi -resolution  California  coastal 
terrain  fields,  respectively. 
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ratio  was  interpolated  exponentially  in  pressure.   U  and  V 
velocity  components  were  interpolated  linearly  in  pressure. 

Topography  was  derived  from  the  U.S.  Navy  global  10 -minute 
elevation  data.  Five  point  horizontally- averaged  smoothing 
was  applied  to  the  elevation  data.  This  provided  a  near- 
realistic  terrain  field  for  the  model  validation  run.  Figure 
1  shows  the  validation  terrain  field. 
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III.  METEOROLOGICAL  SITUATION  (02-03  MAY  1990) 

A.   SYNOPTIC  SCALE  ANALYSIS 

Figures  4,  6,  8,  and  10  show  12 -hour  2.5°  resolution  FNOC 
NOGAPS  500  mb  height  and  temperature  analyses  for  the  period 
0000  UTC  02  May  1990  to  1200  UTC  03  May  1990.  After  initial 
weakening  during  the  first  12  hours,  a  high  amplitude  pressure 
ridge  at  500  mb  in  the  central  East  Pacific  was  slowly 
building  northeast  through  the  period  (Figure  6) .  Subsiding 
northerly  flow  on  the  east  side  of  the  ridge  created  strong 
low- level  and  surface-based  inversions  along  the  entire 
California  coast.  A  low  pressure  cell  over  western  Arizona 
was  stationary  during  the  first  12  hours  of  the  period  and 
deepened  slightly  (Figure  4)  .  During  the  final  24  hours 
(Figures  8  and  10) ,  however,  this  low  filled  and  eventually 
moved  off  to  the  east  - northeast .  A  thermal  trough  over  the 
Pacific  Northwest  at  0000  UTC  02  May  (Figure  4)  moved  south  to 
near  southern  California  by  1200  UTC  03  May  (Figure  10)  .  This 
troughing  was  indicative  of  a  short  wave,  with  its  associated 
positive  vorticity,  in  the  flow  pattern  aloft. 

Figures  3,  5,  7,  and  9  show  12 -hour  2.5°  resolution  FNOC 
NOGAPS  sea  level  pressure  analyses  for  the  period  0000  UTC  02 
May  1990  to  1200  UTC  03  May  1990.  A  high  pressure  cell, 
initially  located  in  the  central  East  Pacific  (Figure  3) ,  had 
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Figure  3.   0000  UTC  02  May  1990  NOGAPS  Sea 
Level  Pressure  (mb)  Analysis. 


Figure  4.  0000  UTC  02  May  1990  NOGAPS  500 
mb  analysis.  Heights  (m,  solid)  and  tem- 
perature (°C,  dashed)  . 
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Figure  5.   12  00  UTC  02  May  199  0  NOGAPS  Sea 
Level  Pressure  (mb)  Analysis. 


Figure  6.  1200  UTC  02  May  1990  NOGAPS  500 
mb  analysis.  Heights  (m,  solid)  and  tem- 
perature (°C,  dashed)  . 
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Figure  7.   0000  UTC  03  May  1990  NOGAPS  Sea 
Level  Pressure  (mb)  Analysis. 


Figure  8.  0000  UTC  03  May  1990  NOGAPS  500 
mb  analysis.  Heights  (m,  solid)  and  tem- 
perature (°C,  dashed)  . 
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Figure  9.   1200  UTC  03  May  1990  NOGAPS  Sea 
Level  Pressure  (mb)  Analysis. 


Figure  10.  1200  UTC  03  May  1990  NOGAPS  500 
mb  analysis.  Heights  (m,  solid)  and  tem- 
perature (°C,  dashed)  . 
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intensified  and  moved  to  near  the  Oregon  coast  by  1200  UTC  03 
May  (Figure  9) .  A  low  pressure  system,  situated  in  southern 
Arizona,  had  moved  slowly  north  and  filled  through  out  the 
period.  Thermal  troughing  in  central  California  was  initially 
coupled  to  this  surface  low,  but  by  1200  UTC  03  May  (Figure  9) 
had  moved  offshore. 

B.   MESOSCALE  ANALYSES 

Several  mesoscale  atmospheric  phenomena  have  been  observed 
and  studied  during  the  period,  0000  UTC  02  May  -1200  UTC  03 
May.  All  three  were  coastal  and  possibly  topographically 
induced  and  are  therefore  relevant  to  this  paper. 

Streed  (1990)  successfully  observed  land/sea  breeze  struc- 
ture for  the  period  02-03  May  1990  using  the  UHF  doppler  wind 
profiler  at  the  Naval  Postgraduate  School,  Monterery,  CA. 
Although  it  is  impossible  to  determine  the  horizontal  extent 
of  the  land- sea  breeze  from  a  single  observation  point, 
vertical  structure  and  intensity  were  determined.  Specifical- 
ly, the  land- sea  breeze  was  confined  to  the  lowest  2500  meters 
of  the  atmosphere.  Although  the  return  limb  of  the  flow  aloft 
is  not  well  defined  for  either  the  land  or  sea  breeze  compo- 
nents, the  lower  limb  is  easily  discerned  in  the  data  (Figure 
11)  to  be  confined  to  the  lowest  1000  meters.  From  0400  UTC 
02  May  1990  to  2000  UTC  02  May  1990,  the  land/valley  breeze 
sets  up  from  the  south,  aligned  with  the  topography  of  the 
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Salinas  River  valley.  The  sea  breeze  persisted  from  approxi- 
mately 2100  UTC  02  May  1990  to  0300  UTC  03  May  1990. 

In  his  case  study  of  a  stratus/fog  event,  Corkill  (1991) 
describes  mesoscale  features  affecting  the  California  region 
during  the  period  0000  UTC  02  May  1990  to  1200  UTC  03  May 
1990.  Surface  mesoscale  analyses  were  made  using  a  multi- 
quadric  interpolation  scheme  (Hardy  199  0) .  The  sea  level 
pressure  analysis  for  0000  UTC  02  May  1990  (Figure  12)  shows 
a  complex  low  pressure  pattern  in  central  California  associat- 
ed with  the  thermal  trough.  Offshore,  flow  is  from  the  north 
with  no  mesoscale  structure  evident.  By  1200  UTC  03  May,  the 
sea  level  pressure  analysis  (Figure  13)  shows  a  mesoscale 
trough  had  developed  offshore  central  California  and  was 
responsible  for  southerly  flow  along  the  coast.  This  situa- 
tion is  similar  to  those  depicted  in  the  Mass  and  Albright 

(1987)  case  studies  of  topographically  trapped  alongshore 
surge  along  the  west  coast  of  North  America.  This  analysis  is 
very  different  from  the  NOGAPS  analysis  (Figure  9) . 

A  third  mesoscale  feature,  the  Catalina  Eddy,  is  discern - 
able  in  the  surface  and  upper  air  data  (see  Chapter  5) .  In 
particular,  soundings  at  San  Clemente  (NSI) ,  El  Toro  (LIO) , 
and  Point  Mugu  (NTD)  suggest  a  closed  Eddy  vortex  at  both  1200 
UTC  02  May  1990  and  at  1200  UTC  03  May  1990.  This  feature  is 
coupled  to  the  southerly  surge  phenomena.   Mass  and  Albright 

(19  89)  describe  the  synoptic  pattern  favorable  for  formation 
of  a  Catalina  Eddy.   The  Catalina  Eddy  results  from  the 
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Figure  12.  0000  UTC  02  May  1990 
Mesoscale  SLP  (mb)  Anal  (Corkill, 
1991)  . 


Figure  13.  1200  UTC  03  May  1990 
Mesoscale  SLP  (mb)  Anal  (Corkill, 
1991)  . 
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interaction  of  synoptic  scale  flow  with  the  formidable 
topographic  barriers  of  the  region.  The  first  process  in  the 
development  of  the  Eddy  is  the  passage  of  a  short-wave  trough 
into  the  Pacific  Northwest  (Figure  6) .  A  pre-existing  850  mb 
trough  in  the  western  United  States  extends  southwestward 
towards  the  southern  California  bight,  while  at  the  surface 
the  interior  heat  trough  extends  northwestward  into  central 
California  (Figure  5) .  These  synoptic  features  intensify  the 
lower  tropospheric  pressure  gradient  along  the  central 
California  coast,  and  strengthen  the  northerly  flow  in  the 
region. 

With  a  strengthening  northerly  flow  approaching  the  San 
Ynez  Mountains  to  the  north  of  the  bight,  synoptic  and 
mesoscale  lee  troughing  become  evident  within  the  bight.  This 
troughing  produces  an  alongshore  pressure  gradient  with  lower 
pressure  to  the  north.  As  a  result,  ageostrophic  southerly 
flow  develops  within  approximately  a  Rossby  radius  (approxi- 
mately 100  km)  of  the  coastal  mountain  barrier.  With  rela- 
tively geostrophic  northerlies  remaining  further  offshore, 
considerable  cyclonic  vorticity  is  created  in  the  coastal 
zone. 

As  the  southerly  flow  increases,  there  is  an  increase  in 
the  depth  of  the  cool,  moist  marine  layer  in  the  coastal  zone, 
resulting  in  a  narrow  coastal  pressure  ridge.  As  the  coastal 
pressure  ridge  intensifies  and  extends  northward,  an  isolated 
low  center  offshore  often  becomes  evident. 
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Catalina  Eddy  events  continue  as  long  as  the  synoptic- 
scale  conditions  support  the  required  alongshore  pressure 
gradient.  As  the  synoptic  pattern  evolves  and  the  inland 
troughing  attenuates,  the  northerly  flow  to  the  west  and  north 
of  Point  Conception  weakens,  the  alongshore  pressure  gradient 
reverses,  and  the  normal  winds  in  the  bight  are  reestablished. 
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IV.  MODEL  PERFORMANCE  FOR  SYNOPTIC  SCALE  PHENOMENA 

A.   GENERAL 

The  model  was  evaluated  for  36  hours  of  integration  during 
the  period  0000  UTC  02  May  to  1200  UTC  03  May.  This  particular 
time  period  was  chosen  for  a  variety  of  reasons.  Not  only  was 
the  Naval  Postgraduate  School  (NPS)  wind  profiler  operating 
during  this  period,  but  the  research  vessel  PT  SUR  was 
offshore  in  Monterey  Bay  providing  upper  air  soundings  every 
six  hours  commencing  at  1800  UTC  02  May.  Thus,  in  addition  to 
routine  NWS  surface  and  upper  air  reports,  two  other  vertical 
high  resolution  data  sets  were  available. 

The  synoptic  situation  was  also  considered  when  choosing 
the  integration  period.  The  synoptic  pattern  was  variable 
enough  during  the  period  to  provide  for  both  alongshore  and 
low  level  off  land  flow  over  the  central  California  coast. 
Additionally,  it  was  responsible  for  generating  significant 
topographic  and  coastal  induced  mesoscale  systems  (see 
previous  chapter) . 

The  evaluation  was  concerned  with  the  model's  performance 
in  predicting  primarily  geopotential  height,  temperature,  and 
wind  fields  and,  secondarily,  moisture  fields.  Model  output 
at  12,  24,  and  36  hour  integrations  were  examined.  Three 
comparison  schemes  were  utilized.   First,  model  output  was 
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compared  with  the  large  scale  interpolated  NOGAPS  analyses 
(with  corresponding  station  plots)  at  the  surface  and  850  mb, 
700  mb,  and  500  mb  constant  pressure  levels.  Second,  cross- 
section  analyses  were  examined.  Third,  vertical  profile 
comparisons  were  made  at  upper  air  sounding  locations. 

Both  model  output  and  interpolated  NOGAPS  fields  were 
plotted  using  the  NPS  Interactive  Digital  Environmental 
Analysis  (IDEA)  Lab  GENAL  general  analysis  package.  IDEA  Lab 
GEMPAK  routines  were  utilized  to  plot  vertical  profiles  and 
perform  model  independent  cross -section  analysis.  The  GEMPAK 
cross -section  analysis  scheme  is  based  upon  cubic  spline 
interpolation  of  sounding  data.  This  analysis  scheme  is 
limited  in  its  ability  to  resolve  small  scale  features,  and 
its  interpolation  scheme  often  generates  fictious  mesoscale 
structure,  particularly  in  the  higher  elevations.  For  these 
reasons,  it's  primary  use  was  to  verify  the  model's  handling 
of  large  scale  processes. 

The  GEMPAK  cross -section  analysis  is  done  on  a  straight 
line  between  the  first  and  last  stations;  the  positions  of 
intervening  stations  are  proportional  to  the  perpendicular 
projections  of  the  actual  positions  onto  the  section  line. 
Since  at  least  four  stations  are  required  for  the  scheme,  four 
of  the  section  lines  were  extended  beyond  the  model  domain. 
Figure  14  shows  the  section  lines;  because  of  data  nonavail- 
ability, not  all  section  lines  were  used  at  each  evaluation 
time.   Orientation  of  the  sections  are  based  on  sounding 
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station  alignment  rather  than  concern  with  examining  a 
specific  atmospheric  phenomena. 

Input  fields  at  0000  UTC  02  May  and  output  fields  at  1200 
UTC  02  May  are  discussed  below  in  detail.  Supporting  charts 
are  included  for  these  times  in  order  to  illustrate  the 
validation  process.  The  same  validation  process  was  applied 
to  the  remaining  two  times,  0000  and  1200  UTC  03  May;  however, 
only  a  summary  and  any  significant  departures  during  these 
periods  from  previously  discussed  model  performance  is 
included. 

B.   ANALYSIS  AT  0000  UTC  02  MAY  1990  (MODEL  INPUT  FIELDS) 

In  order  to  determine  the  fidelity  of  the  initial  large- 
scale  analysis,  comparison  is  made  of  the  0000  UTC  02  May  1990 
interpolated  NOGAPS  fields  (which  serve  as  the  model  initial- 
ization fields)  to  observations  and  GEMPAK  cross-section 
analyses . 

1.   Constant  Pressure  Surfaces 
a .  Surface 

The  vertical  interpolation  scheme  uses  an  algorithm 
which  extrapolates  temperature  from  the  fifth  lowest  sigma 
level  downward  to  obtain  surface  pressure.  This  leads  to 
disparity  between  interpolated  surface  pressure  fields  and  the 
large-scale  NOGAPS  sea  level  pressure  pattern.  In  general, 
however,  the  interpolated  NOGAPS  sea  level  pressure  pattern 
appears  reasonable  (Figure  15) .   A  single  ship  report  (XCM1) 


26 


H    <D 

rt  u 

-H    3 

4J  iJ 

-H    (tJ 

C    ^ 

H    Q) 

& 

o  B 

<J\   Q) 

CT>  H 

H 

T) 

>iC 

rt  (C 

H  _ 

SI 

■^— • 

u 

S2 

3 

O    CD 

O    CQ 

O    CD 

CQ 

O    U 

a* 

-H 
CO 

e   15. 

evel 

>1 

S 

U  J 

-H    CD 

0 

fr.  OD 

27 


suggests  slightly  weaker  pressure  gradients  offshore  than 
analyzed.  The  interpolation  has  produced  slightly  stronger 
ridging  over  the  mountains  than  indicated  on  the  synoptic 
scale  analysis  (Figure  3)  .  Thermal  troughing  over  central 
California  appears  reasonable  in  orientation  and  magnitude. 
West  and  southwest  winds  at  several  coastal  stations  in 
central  and  southern  California  suggest  a  more  complex  mesos- 
cale  structure  exists  than  indicated  by  the  interpolated 
analysis.  Sea  surface  temperatures  look  reasonable  based  on 
climatology  and  generally  reflect  the  expected  pattern  for  a 
cold  southward  flowing  California  current. 
b.  850  Mb  Constant  Pressure 

The  850  mb  height  and  temperature  analysis  (Figure 
16)  does  not  deviate  significantly  from  the  synoptic  scale 
analysis  (not  shown) .  Gaps  in  the  contours  indicate  where  the 
850  mb  pressure  level  intersects  topography;  no  extrapolated 
analysis  below  the  terrain  was  done  within  that  region.  Weak 
thermal  ridging  with  associated  warm  air  advection  exists 
offshore  of  central  California;  observations  support  this 
thermal  analysis.  The  initial  wind  field,  shown  in  Figure  17 
by  a  wind  vector  at  every  other  point  of  the  model  grid, 
departs  significantly  from  observed  winds,  particularly  at 
Winnemuca  (WMC) ,  Medford  (MFR)  and  Desert  Rock  (DRA) .  These 
initial  wind  field  deviations  most  likely  reflect  observed 
topographic  effects. 
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The  mixing  ratio  analysis  at  850  mb  shows  considerable  dry- 
air  off  the  northern  and  central  California  coast.  A  moisture 
maximum  exists  over  western  Nevada  east  of  the  Sierra  Nevada. 
Examination  of  lower  level  moisture  fields  show  this  maximum 
tilts  eastward  and  weakens  with  height.  As  previously  men- 
tioned, there  is  a  significant  departure  between  the  NOGAPS 
and  NMC  GDAS  moisture  analyses.  This  difference  is  most 
notable  at  850  mb  where  NOGAPS  and  GDAS  moisture  analyses 
appear  almost  180  degrees  out  of  phase.  Sounding  data  and 
cross -section  analysis  suggest  that  the  NMC  scheme  may  be  more 
correct  in  handling  moisture. 

c.  700  Mb  Constant  Pressure 

The  700  mb  pressure  surface  lies  just  above  the 
highest  terrain  of  the  Sierra  Nevada.  At  this  level,  the 
atmosphere  is  nearly  geostrophic.  Wind,  height,  and  thermal 
observations  should  therefore  be  consistent  with  the  initial 
analysis  fields.  The  height  and  temperature  analysis  (figures 
not  shown)  shows  only  minimal  noise  over  the  mountains.  This 
noise  is  an  artifact  of  the  vertical  interpolation  scheme  and 
reflects  extrapolation  of  the  temperature  field  below  the 
lowest  sigma  level  to  the  surface.  The  noise  is  most  pro- 
nounced at  low- levels  and  over  mountainous  terrain.  The 
height  field  and  sounding  data  correlate  well  at  all  loca- 
tions. The  thermal  pattern  generally  reflect  the  observa- 
tions.   There  is  some  problem,   however,   in  positioning 
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specific  temperature  contours,  most  notably  in  the  weak 
thermal  gradient  over  Nevada.  Wind  field  and  observations  are 
in  fair  agreement  everywhere  but  at  WMC.  The  moisture 
analysis  shows  the  western  Nevada  maximum  as  well  as  increas- 
ing the  moisture  to  the  northwest.  Drier  air  lies  off  the 
central  and  southern  California  coasts.  Considerably  greater 
agreement  exists  between  the  NOGAPS  and  GDAS  moisture  fields 
at  this  level  than  at  850  mb  and  below. 

d.  500  Mb  Constant  Pressure 

At  500  mb,  terrain  influences  are  non-existent. 
Height,  thermal,  (Figure  18)  and  wind  fields  (Figure  19)  are 
in  good  agreement  with  observed  data.  Because  of  the  low 
moisture  content  in  the  atmosphere  at  these  levels,  the 
moisture  field  was  not  verified. 
2.   CROSS -SECTION  ANALYSES 

Cross -section  analyses  along  the  Medford  (MFR)  to  San 
Nicholas  Island  (NSI)  transect  (see  Figure  14)  were  made  for 
both  the  interpolated  analysis  fields  (using  GENAL)  and  upper 
air  sounding  data  (using  GEMPAK) .  As  expected,  no  BL  struc- 
ture is  evident  in  the  interpolated  fields.  Potential 
temperature  cross-sections  (Figures  20  and  21)  at  0000  UTC  02 
May  1990  agree  on  general  patterns,  although  the  GEMPAK 
analysis  does  suggest  some  BL  structure  south  of  Oakland. 
Considering  the  time  of  day  (late  afternoon) ,  BL  instability 
is  expected. 
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Figure  20.   0000  UTC  02  May  1990  MFR-NSI  Theta 
(°K)  Model  Input  Cross -section  Analysis. 


Figure  21.   0000  UTC  02  May  1990  MFR-VBG  Theta 
(°K)  GEMPAK  Cross -section  Analysis. 
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Mixing  ratio  cross-sections  (Figures  22  and  23)  disagree 
significantly,  particularly  south  of  OAK.  The  GEMPAK  analysis 
shows  considerably  drier  air  above  900  mb,  suggesting  that  the 
model  input  data  may  be  too  moist  in  this  region.  The  NMC 
GDAS  moisture  analysis  lends  support  to  this  conclusion. 

Both  the  u-velocity  cross-sections  (Figures  24  and  25)  and 
v-velocity  cross-sections  (Figures  26  and  27)  above  approxi- 
mately 800  mb  agree  well,  but  show  some  disparity  at  lower 
levels.  Most  notably,  the  lower  GEMPAK  v- component  velocities 
in  the  vicinity  of  MFR  reflect  topographic  influences  not 
discerned  in  the  interpolated  field  analysis.  Also,  a 
southerly  wind  component  near  Vandenberg  AFB  (VBG)  is  not 
analyzed  in  the  interpolated  fields. 

C.   MODEL  OUTPUT  AT  12  00Z  02  MAY  1990 

In  order  to  determine  the  model's  ability  to  successfully 
forecast  large-scale  synoptic  features,  comparison  is  now  made 
of  the  12  hour  model  simulation  at  1200  UTC  02  May  1990  to  the 
corresponding  interpolated  analyses,  observations,  and  GEMPAK 
cross -section  analyses. 

1.   Constant  Pressure  Surfaces 
a .  Surface 

The  model's  forecast  surface  pressure  and  tempera- 
ture fields  (Figure  28)  differs  significantly  from  the 
interpolated  analysis  (Figure  29)  .  Most  notably,  the  model 
depicts  a  large  region  of  broad  scale  lower  pressure  off  the 
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Figure  22.   0000  UTC  02  May  1990  MFR-NSI  Mixing 
Ratio  (g/kg)  Model  Input  Cross -section  Analysis. 


Figure  23.   0000  UTC  02  May  1990  MFR-VBG  Mixing 
Ratio  (g/kg)  Model  Input  Cross-section  Analysis . 
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Figure  24.  0000  UTC  02  May  1990  MFR-NSI  U- 
velocity  Model  Input  Cross-section  Analysis, 
Eastward  (m/s,  solid).   Westward  (m/s,  dashed). 


Figure  25.  0000  UTC  02  May  1990  MFR-VBG  U- 
velocity  Model  Input  Cross- section  Analysis. 
Eastward  (m/s,  solid) .   Westward  (m/s,  dashed) . 
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Figure  26.  0000  UTC  02  May  1990  MFR-NSI  V- 
velocity  Model  Input  Cross-section  Analysis, 
Northward  (m/s,  solid) .  Southward  (m/s,  dashed) , 


500 


J  000' 


Figure  27.  0000  UTC  02  May  1990  MFR-VBG  V- 
velocity  GEMPAK  Cross -section  Analysis.  North- 
ward (m/s,  positive) .  Southward  (m/s,  negative) , 
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southern  California  coast  whereas  the  analysis  shows  weak 
ridging  through  the  same  area.  A  lack  of  data  in  the  region 
precludes  verifying  either  the  model  or  analysis. 

Over  land,  as  expected,  the  model  shows  considerable  more 
structure.  Low  pressure  regions  over  the  Santa  Ynez  and 
Sierra  Nevada  Mountains,  however,  appear  to  be  artifacts  of 
the  model's  sea  level  pressure  conversion  algorithm  and  are 
not  supported  by  station  data.  Within  the  San  Joaquin  Valley, 
pressure  gradient  is  too  flat  to  discern  verifiable  features. 
Elsewhere,  model  output  appears  reasonable. 

Temperature  fields  over  water  are  nearly  identical  for 
both  the  model  and  interpolated  analysis.  Over  land,  model 
surface  temperatures  also  look  reasonable.  As  expected  for 
12  00  UTC  early  morning  hours,  ground  temperature  is  signifi- 
cantly cooler  than  SST.  The  highest  surface  temperatures  are 
found  at  low  elevations  (in  the  valleys) ,  lowest  temperatures 
at  the  highest  elevations. 

Jb.  850  Mb  Constant  Pressure 

Model  height  fields  (Figure  30)  show  the  same  broad 
scale  troughing  over  the  south  central  portion  of  the  domain 
as  the  surface  pressure  output.  Observed  heights  (Figure  31) 
at  VBG,  NTD,  LIO  and  NSI  suggest  the  interpolated  analysis 
more  closely  represents  actual  height  fields,  as  expected  from 
a  00 -hour  analysis  versus  a  12 -hour  forecast.   The  model's 
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temperature  field,  however,  appears  to  more  closely  represent 
observed  conditions  than  the  interpolated  analysis. 

Wind  fields  (Figures  32  and  33)  are  better  represented  by 
model  output  over  rough  terrain  (NID,  DRA,  WMC)  whereas  the 
analysis  is  better  along  the  coast  (OAK,  VBG,  NTD,  LIO,  and 
NSI)  . 

c.  700  Mb  Constant  Pressure 

In  general,  model  and  analysis  height  fields 
(Figures  34  and  35)  agree  both  with  each  other  and  with 
observations.  The  only  significant  difference  lies  in  the 
southern  part  of  the  domain  over  water  where  the  analysis 
shows  more  ridging  than  the  model.  Again,  little  data  is 
available  to  support  either  representation.  With  regard  to 
the  temperature  fields,  neither  the  analysis  nor  model  can  be 
preferred  based  solely  on  station  data.  However,  as  at  850 
mb,  the  model  shows  much  more  temperature  structure  than  the 
analysis . 

Winds  fields  (Figures  36  and  37)  show  two  major  differenc- 
es. The  analysis  shows  north- northwest  flow  in  the  Califor- 
nia bight,  the  model  northerly  flow.  This  disparity  reflects 
the  difference  in  height  fields  between  the  two  grids.  The 
second  difference  lies  inland  over  the  Sierra  Nevada  where  the 
model  shows  topographically  induced  convergence.  Again,  the 
lack  of  data  precludes  a  determination  of  whether  the  model  or 
analysis  best  represent  actual  conditions. 
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d.  500   Mb  Constant  Pressure 

The  problem  of  height  field  (Figures  3  8  and  39) 
alignment  in  the  southern  part  of  the  domain  still  persists, 
but  to  a  lesser  extent.  Both  model  and  analysis  show  good 
agreement  in  temperature  (Figures  38  and  39)  and  wind  fields 
(Figures  40  and  41).  As  expected  when  comparing  12-hour 
forecast  fields  to  analysis  fields,  the  observed  winds  at  VBG, 
NSI,  and  NID  suggest  that  the  analysis,  rather  than  the  model, 
better  represents  actual  conditions. 

2.   Cross -section  Analyses 
a.    MFR    to  NSI 

Figures  42  and  43  show  GEMPAK  and  model  potential 
temperature  cross -sections,  respectively.  The  most  important 
feature  in  the  GEMPAK  analysis  is  a  strong  stable  boundary 
layer  below  a  thermal  inversion  near  900  mb  from  north  of  OAK 
to  NSI.  This  feature  is  successfully  replicated  by  the  model. 

There  is  also  a  general  similarity  in  the  wind  patterns 
depicted  by  both  the  GEMPAK  and  model  cross  -  sections .  With 
regard  to  the  u- component  of  the  wind  (Figures  44  and  45)  , 
both  cross-sections  show  westerlies  below  the  inversion. 
Above  the  inversion,  the  model  successfully  forecasts  the 
easterly  wind  maximum  over  OAK  near  9  00  mb.  Further  south, 
towards  VBG  and  NSI,  the  analysis  shows  a  westerly  component 
aloft  between  700  and  800  mb  whereas  the  model  only  hints  at 
this  feature. 
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Figure  42.   1200  UTC  02  May  1990  MFR-NSI  Theta 
(°K)  Model  Output  Cross -section  Analysis. 


Figure  43.   1200  UTC  02  May  1990  MFR-NSI  Theta 
(°K)  GEMPAK  Cross -section  Analysis. 
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Figure  44.  1200  UTC  02  May  1990  MFR-NSI  U- 
velocity  Model  Output  Cross -section  Analysis, 
Eastward  (m/s,  solid).   Westward  (m/s, dashed) . 


Figure  45.    1200  UTC  02  May  1990  MFR-NSI  U- 
velocity  GEMPAK  Cross -section  Analysis.  Eastward 
(m/s,  positive).   Westward  (m/s, negative) . 
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The  v- component  cross -sections  (Figures  46  and  47)  also 
are  generally  similar.  Northerly  flow  is  prevalent  at  all 
levels.  Magnitudes  are  reasonably  close;  however,  the  model 
predicts  two  wind  maxima  over  OAK,  one  below  950  mb  in  the 
stable  layer,  the  other  aloft  near  700  mb.  The  GEMPAK 
analysis  shows  only  one  wind  maximum  near  875  mb. 

Both  model   (Figures  48)  and  GEMPAK  moisture  analysis 

(Figure  49)  agree  fairly  well  both  in  pattern  and  magnitude. In 

particular,  both  depict  more  moisture  at  low- levels;   aloft, 

both  show  the  general  trend  of  moisture  decreasing  north  to 

south. 

b.    NSI    to  ELY 

Dramatic  changes  in  topography  along  this  cross - 
section  preclude  a  complete  GEMPAK  analysis  at  lower  levels. 
Still,  enough  of  the  analysis  is  complete  to  allow  evaluation 
of  model  performance. 

Model  (Figure  50)  and  GEMPAK  (Figure  51)  potential 
temperature  analyses  show  comparable  structure  above  the  PBL. 
Within  the  PBL,  the  GEMPAK  analysis  hints  at  a  strong  stable 
layer  west  of  the  mountains  between  NSI  and  EDW  and  over  the 
high  plains  between  DRA  and  ELY.  This  layer  is  replicated  by 
the  model . 

U- component  cross -sections  (Figures  52  and  53)  generally 
agree  above  800  mb  and  to  the  east  at  lower  levels,  but  three 
differences  are  worthy  of  note.   First,  the  model  does  not 
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Figure  46.  1200  UTC  02  May  1990  MFR-NSI  V- 
velocity  Model  Input  Cross-section  Analysis, 
Northward  (m/s, solid) .   Southward  (m/s,  dashed). 


Figure  47.  1200  UTC  02  May  1990  MFR-VBG  V- 
velocity  GEMPAK  Cross -section  Analysis.  North- 
ward (m/s,  positive) .  Southward  (m/s,  negative) . 


51 


Figure  48.   1200  UTC  02  May  1990  MFR-NSI  Mixing 
Ratio  (g/Jcg)  Model  Output  Cross -section  Analysis . 


Figure  49.   1200  UTC  02  May  1990  MFR-NSI  Mixing 
Ratio  (g/kg)  GEMPAK  Cross-section  Analysis. 
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Figure  50.    1200  UTC  02  May  1990  NSI-ELY  Theta 
(°K)  Model  Output  Cross -section  Analysis. 
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Figure  51.    1200  UTC  02  May  1990  NSI-ELY  Theta 
(°K)  GEMPAK  Cross -section  Analysis. 
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Figure  52.  1200  UTC  02  May  1990  NSI-ELY  U- 
velocity  Model  Output  Cross -section  Analysis. 
Eastward  (m/s,  solid) .   Westward  (m/s,  dashed) . 


Figure  53.    1200  UTC  02  May  1990  NSI-ELY  U- 
velocity  GEMPAK  Cross -sect  ion  Analysis.  Eastward 
(m/s,  positive) .   Westward  (m/s,  negative) . 
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fully  estimate  the  westerly  component  at  NSI  below  800  mb.  It 
shows  westerlies  only  at  very  low  levels  (below  950  mb)  , 
although  it  tries  to  capture  the  pattern  by  depicting  light 
easterlies  in  that  region.  Second,  the  GEMPAK  analysis  shows 
a  westerly  maximum  over  NID  at  800  mb.  This  wind  maximum  is 
the  result  of  a  misreported  wind  direction  at  that  level. 
Last,  the  model  shows  a  region  of  strong  low  level  easterlies 
along  the  windward  slope  of  the  coastal  mountain  range.  No 
data  is  available  to  verify  the  existence  of  this  feature. 

V- component  cross -sections  (Figures  54  and  55)  are  in  good 
agreement.  The  general  pattern  of  lower  wind  speeds  to  the 
east  and  west  is  adequately  depicted  by  the  model.  No  major 
disparities  are  evident. 

The  model's  moisture  cross -section  shows  considerably 
higher  mixing  ratio  values  than  analyzed,  possibly  a  result  of 
poor  initial  moisture  input  fields. 

3.   Vertical  Profiles 

Vertical  profiles  from  all  upper  air  observations  made 
within  the  model's  domain  at  1200  UTC  02  May  1990  were 
compared  with  vertical  profiles  obtained  from  model  output  at 
their  locations.  Figures  56  through  59  show  observed  vertical 
profiles  versus  profiles  of  model  output  for  two  of  these 
upper  air  reporting  locations,  NSI  and  ELY.  These  sites  not 
only  display  typical  model  output  characteristics  but  are 
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Figure  54.  1200  UTC  02  May  1990  NSI-ELY  V- 
velocity  Model  Output  Cross -section  Analysis. 
Northward  (m/s,  solid) .  Southward  (m/s,  dashed) . 


Figure  55.  1200  UTC  02  May  1990  NSI-ELY  V- 
velocity  GEMPAK  Cross -section  Analysis.  North- 
ward (m/s,  positive) .  Southward  (m/s,  negative) . 
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Figure  56.   1200  UTC  02  May  1990 
NSI  Observed  Vertical  Profile. 
(Pressure,  mb)  .  (Temperature  and 
Dew  Point,  °C) 
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Figure  57.   1200  UTC  02  May  1990 
Model  Vertical  Profile  at  NSI. 
(Pressure,  mb)  .  (Temperature  and 
Dew  Point,  °C)  . 
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Figure  58.   1200  UTC  02  May  1990 
ELY  Observed  Vertical  Profile. 
(Pressure,  mb)  .   (Temperature  and 
Dew  Point,  °C)  . 


Figure  59.   1200  UTC  02  May  1990 
Model  Vertical  Profile  at  ELY. 
(Pressure,  mb)  .   (Temperature  and 
Dew  Point,  °C)  . 
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representative  of  the  synoptic  and  geographic  diversity 
within  the  model's  domain. 

The  observed  temperature  profile  at  NSI  shows,  as  expect- 
ed, a  strong  elevated  subsidence  inversion  with  a  shallow 
stratus  deck  below.  The  model  misses  this  feature  and  indi- 
cates a  comparable  strength  surface-based  inversion  and  no 
stratus  deck.  Above  the  inversion,  the  model  replicates  well 
the  vertical  thermal  structure.  Moisture  profiles  are  handled 
poorly  by  the  model.  The  model  fails  to  dry  out  the  atmo- 
sphere at  the  top  of  the  inversion  or  to  depict  finer  scale 
features  aloft.  Except  for  the  winds  directly  above  the 
inversion,  model  wind  profiles  are  within  5  kts  and  20  degrees 
of  observed  winds  at  all  levels. 

The  model  is  also  evaluated  at  ELY,  an  inland  mountain 
station.  As  with  the  NSI,  temperature  profiles  are  adequately 
handled  by  the  model,  although  the  model  misses  completely  the 
observed  nocturnal  surface-based  inversion.  Moisture  profiles 
are  marginal  and  correlate  only  on  the  largest  scale.  Model 
and  observed  wind  profiles  differ  by  less  than  20  degrees  and 
5  knots  at  all  levels. 

D.   MODEL  OUTPUT  AT  0000  UTC  03  MAY  1990 

At  0000  UTC  03  May  1990,  low-level  model  (surface  and  850 
mb)  pressure/height  fields  show  discrepancies  with  the 
analyses  on  their  orientation  in  the  California  bight  region. 
The  model  overf orecasts  the  intensity  of  the  low  pressure 
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system  which  dominates  most  of  coastal  southern  California. 
Unlike  at  1200  UTC  02  May  1990,  significant  surface  pressure 
differences  between  model  output  and  both  observations  and 
analyzed  fields  (minus  4-5  mb)  are  evident  within  the  bight 
region.  Though  not  as  dramatic,  forecasted  height  values  at 
850  mb  differ  by  minus  15-20  meters  from  analyzed  values. 

Forecast  temperature  fields  remain  the  model's  strong 
suite.  Very  little  difference  was  noted  between  model  output 
and  observed  temperatures,  regardless  of  location  or  height. 
This  is  particularly  noteworthy  considering  the  diurnal 
heating  effects  which  have  occurred  during  the  12  hour 
integration  period. 

The  model  also  performs  well  with  regard  to  wind  fields  at 
all  levels.  Although  only  eight  vertical  soundings  were 
available  to  verify  model  output,  no  significant  difference 
between  model  and  observed  winds  were  found. 

The  model  continues  to  forecast  an  overabundance  of 
moisture  at  all  levels  in  the  atmosphere.  As  at  1200  UTC  02 
May,  the  spatial  distribution  pattern  is  similar  to  the 
analyses,  however,  the  quantity  of  moisture  predicted  far 
exceeds  the  analyzed  values  at  nearly  all  locations. 

E.   MODEL  OUTPUT  AT  1200  UTC  03  MAY  1990 

At  the  surface,  the  model  continues  to  overdeepen  the  low 
pressure  system  in  central  and  southern  California.  Differ- 
ences of  minus  6-7  mbs  are  noted  along  the  central  California 
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coast,  although  the  contour  pattern  and  alignment  in  that 
region  correlates  well  with  the  analysis.  The  problem  with 
contour  alignment  in  the  southern  latitudes  of  the  domain, 
however,  persists. 

At  850  mb,  observed  winds  and  forecast  wind  fields  are  in 
fair  agreement,  with  closer  agreement  over  the  less  mountain- 
ous regions.  Geopotential  height  fields,  like  that  of  surface 
pressure,  are  significantly  lower  than  the  analysis  (15-20 
meters) .  Above  850  mb,  the  difference  between  forecast  height 
fields  and  analyzed  heights  diminish  with  elevation. 

Again,  temperature  fields  are  well  forecasted  by  the 
model.  Both  vertical  cross-sections  and  station  profiles  show 
fairly  good  agreement  at  all  levels,  although  the  model  was 
not  able  to  fully  capture  inversions  as  well  as  it  had  for  the 
previous  24 -hour  integration  time. 

No  change  in  the  model's  general  handling  of  moisture  was 
seen  at  1200  UTC  03  May. 


61 


V.   MODEL  PERFORMANCE  FOR  MESOSCALE  PHENOMENA 

A.   GENERAL 

The  purpose  of  the  model  assessment  is  to  evaluate  the 
model's  potential  for  use  in  studies  of  mesoscale  coastal  and 
topographically  influenced  meteorological  features.  Following 
the  assessment  of  the  model's  capability  to  replicate  synoptic 
scale  systems  in  the  last  chapter,  an  evaluation  was  made  of 
the  model's  capability  to  forecast  observed  mesoscale  flow. 

Three  types  of  mesoscale  atmospheric  phenomena  were 
observed  along  the  California  coast  during  the  period  0000  UTC 
02  May  1990  -  1200  UTC  03  May  1990.  These  phenomena  include 
sea  and  land  breezes  in  the  vicinity  of  Monterey,  topographi- 
cally trapped  density  currents  (the  southerly  surge) ,  and  the 
Catalina  Eddy.  As  discussed  in  Chapter  3,  all  three  are 
coastal  phenomena  which  are  heavily  influenced  by  local 
topography.  Analyses  of  each  phenomena  were  made  using 
available  observations;  however,  these  observations  were 
limited  and  in  each  case  the  mesoscale  features  could  not 
fully  be  described.  Still,  the  analyses  were  sufficiently 
detailed  to  allow  viable  comparison  with  model  output  and  thus 
provide  a  means  to  assess  model  mesoscale  performance. 
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B.   LAND/SEA  BREEZES  AT  PORT  ORD#  MONTEREY,  CALIFORNIA 

Figure  60  is  a  time-height  series  for  the  observed  (dashed 
lines)  and  model  output  (solid  lines)  of  u-wind  components  at 
the  NPS  doppler  wind  profiler  site  at  Fort  Ord,  California, 
during  the  period  0000  UTC  02  May  1990  -  1200  UTC  03  May  1990. 
Although  the  observed  winds  are  a  combination  of  flows  at  all 
scales  (i.e.,  synoptic,  mesoscale,  local),  information  about 
the  lower  limb  of  the  sea  breeze  can  be  discerned  below  1400 
meters.  In  particular,  observed  westerly  wind  maxima  (8-10 
m/s)  at  0000  UTC  02  May  and  0000  UTC  03  May  (1600  LST) 
indicate  the  sea  breeze  component  of  the  circulation.  The 
upper  limb  is  masked  by  larger  scale  flow  and  not  discernable. 

Model  output  (solid  lines)  does  quite  well  in  replicating 
the  time  of  the  sea  breeze  component  of  the  circulation.  Both 
observed  0000  UTC  02  May  1990  and  0000  UTC  03  May  westerly 
wind  maxima  are  captured  by  the  model,  although  the  magnitude 
of  the  model  output  is  considerably  less  (2-4  m/s)  than  the 
observed  values.  This  difference,  however,  may  be  caused  by 
small  deviations  between  observations  and  model  output  in  the 
larger  scale  flow  field. 

Above  14  00  meters,  the  model  shows  some  disagreement  with 
observed  winds.  These  disparities  are  largely  due  by  small 
differences  (10-15  degrees)  between  model  wind  direction  and 
observed  wind  direction  aloft. 
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The  land  breeze  at  the  profiler  site  sets  up  with  the 
alignment  of  the  Salinas  River  Valley  from  the  southeast. 
The  observed  v- component  time-height  series  (Figure  61,  dashed 
lines)  clearly  shows  southerly  wind  maxima  near  the  ground  at 
1200  UTC  (0400  LST)  02  May  and  0900  UTC  (0100  LST)  03  May. 
These  maxima  probably  indicate  the  combined  effects  of  both 
the  land  breeze  and  the  southerly  surge,  with  the  strongest 
southerly  surge  component  most  likely  occurring  at  0900  UTC  03 
May.  It  is  impossible  to  clearly  identify  only  the  land 
breeze  component  of  the  data  so  no  conclusive  verification  of 
the  model's  land  breeze  can  be  made.  However,  the  strong 
correlation  of  the  timing  in  diurnal  variations  between 
observation  and  model  output  (solid  lines)  clearly  shows  the 
model  has  captured  some  of  the  land  breeze  characteristics. 

C.   THE  SOUTHERLY  SURGE  ALONG  THE  CALIFORNIA  COAST 

Surface  southerlies,  associated  with  a  weak  southerly 
surge,  occurred  during  the  period  0000  UTC  02  May  -  1200  UTC 
03  May  1990.  These  southerlies  are  diurnal  in  nature,  with 
the  strongest  observed  flow  occurring  at  1200  UTC  (0400  LST) 
daily.  (Modle  output  indicates  strongest  southerly  coastal 
flow  occurs  at  1800  UTC  02  May  1990.  No  data  is  available, 
however,  to  confirm  this  output.)  This  diurnal  pattern 
suggests  a  link  between  the  strength  of  the  southerly  coastal 
flow,  the  Catalina  Eddy  (which  also  exhibits  similar  diurnal 
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influences) ,  the  land/sea  breeze  along  the  coast,  and  possibly 
the  strength  of  the  San  Joaquin  Valley  afternoon  thermal  low. 

Corkill  (1991) ,  utilizing  a  mesoscale  analysis  scheme 
based  on  multi-quadric  interpolation,  shows  the  strongest 
surface  southerlies  occurred  at  1200  UTC  03  May.  Figure  3  is 
his  sea  level  pressure  analysis  for  that  time.  A  narrow  band 
of  high  pressure  extends  along  the  California  coast  from 
offshore  near  Point  Conception  to  Monterey  Bay.  As  indicated 
in  Chapter  3,  southerly  flow  leads  to  increased  marine  layer 
depth.  Damming  of  this  cool  marine  air  by  the  coastal 
mountain  range  leads  to  formation  of  this  mesoscale  high 
pressure  feature.  Vertical  wind  profiles  at  Vandenberg  (VBG) , 
from  RV  PT  SUR,  and  at  the  NPS  profiler  site  provide  addition- 
al information  on  the  structure  of  the  surge  phenomena  at 
higher  levels. 

Figure  62  is  the  model's  sea  level  pressure  forecast  for 
1200  UTC  03  May  1990.  The  model  indicates  weak  coastal 
pressure  ridging  just  offshore  between  Point  Conception  and 
Monterey  Bay.  The  model's  southerly  flow  is  therefore 
considerably  weaker  than  the  observed  flow.  There  are  several 
explanations  for  this  discrepancy  in  the  model's  output. 

First,  note  that  the  model  surface  pressures  are  consider- 
ably deeper  than  observations,  an  artifact  of  the  model's 
incorrect  forecast  of  the  synoptic  height  fields  discussed  in 
the  previous  chapter.  This  incorrect  forecast  precludes  good 
model  replication  of  the  coastal  southerlies.  As  discussed  in 
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Chapter  3,  an  alongshore  pressure  gradient  with  lower  pressure 
to  the  north  is  required  for  ageostrophic  coastal  southerly- 
flow  development.  The  model's  overdeepening  of  low  pressure 
to  the  south  creates  only  a  weak  northward  alongshore  pressure 
gradient.  Hence,  the  model's  southerly  flow  is  not  developed 
to  the  same  extent  as  the  observed  southerly  flow. 

Second,  the  modeled  coastal  terrain  is  considerably 
smoother  than  the  actual  terrain  features.  Gaps  in  the  model 
topography  may  have  "leaked"  energy  through  the  coastal  range 
into  the  interior  San  Joaquin  Valley  region.  The  modeled  cold 
air  damming  would  not  be  as  strong,  therefore,  as  what 
actually  occurred,  and  the  modeled  southerly  flow  would  be 
considerably  weaker  than  observed. 

Finally,  the  inability  of  the  model  to  adequately  develop 
the  marine  layer  may  have  contributed  significantly  to  the 
model's  underf orecast  of  the  southerly  surge's  intensity. 

Interestingly,  the  observed  VBG  sounding  shows  no  souther- 
ly wind  components  at  any  level  during  the  entire  period  0000 
UTC,  02  May- 12  00  UTC,  03  May.  This  suggests  that  the  ridge 
lies  west  of  VBG  throughout  the  entire  period. 

D.   TOPOGRAPHICALLY  INDUCED  GRAVITY  WAVES 

Observed  winds  aloft  at  the  NPS  profiler  site  (Figure  61) 
show  relatively  strong  southerlies  at  0600  UTC  02  May  and  06 
UTC  03  May  from  the  surface  to  nearly  3500  meters  (700  mb) . 
The  flow  below  1000  meters  is  most  likely  associated  with  the 
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land  breeze.  There  is  a  fairly  good  correlation  between  the 
southerly  wind  maxima  near  3000  meters  at  0600  UTC  and  the 
model's  replication  of  the  flow  trend  (weakened  northerly 
flow)  at  those  same  levels  and  times. 

Although  no  definitive  conclusions  can  be  made,  the  height 
of  these  observed  southerly  wind  maxima  most  likely  precludes 
their  association  with  the  low- level  southerly  flow  field. 
Both  OAK  0000  UTC  03  May  and  RV  PT  SUR  0600  UTC  03  May  1990 
soundings  show  a  weak  inversion  and  significant  drying  out  at 
that  height.  These  observations,  along  with  the  diurnal 
recurrence  of  the  southerly  maxima,  indicate  that  they  might 
be  gravity  waves  on  an  internal  boundary  layer  at  that  level. 
This  internal  boundary  layer  could  be  generated  by  broad  scale 
northeasterly  flow  off  the  Sierra  Nevada  which  manifests 
itself  on  the  west  coast  after  the  thermal  low  in  the  San 
Joaquin  Valley  is  sufficiently  weak  to  allow  energy  to 
propagate  westward.  A  similar  explanation  was  offered  by 
Streed  (1990)  who  suggested  these  anomalies  were  generated  by 
lower,  local  topographic  features. 

E.   THE  CATALINA  EDDY 

Mass  and  Albright  (19  89)  describe  the  full  evolution  of  a 
composite  Catalina  Eddy  event.  The  Catalina  Eddy  results  from 
the  interaction  between  synoptic  scale  flow,  regional  diurnal 
circulations,  and  the  complex  topography  surrounding  the 
California  bight.    Therefore  the  model's  handling  of  the 
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Catalina  Eddy  which  occurred  during  the  period  0000  UTC  02  May 
1990  -  1200  UTC  03  May  1990  is  of  particular  interest  when 
assessing  the  model's  performance. 

Aside  from  surface  observations,  only  soundings  from  LIO, 
NSI,  VBG,  and  NTD  were  available  to  define  the  Eddy  structure. 
These  profiles  were  linearly  interpolated  to  obtain  values  at 
model  levels.  The  model's  overprediction  of  low  pressure  in 
the  bight  has  already  been  addressed;  therefore,  only  the 
model's  Eddy  wind  fields  will  be  evaluated. 

At  12  00  UTC  02  May  199  0,  soundings  at  NSI,  LIO  and  NTD 
indicate  the  Eddy  vortex  was  located  near  Santa  Catalina 
Island.  Based  on  these  observations,  a  closed  vortex  was 
suggested  at  all  heights  up  to  920  mb.  The  model  replicates 
the  depth  of  the  vortex,  although  it  places  the  vortex  center 
farther  south  than  the  analyzed  position.  Figure  63  shows 
model  output  wind  fields  at  9  80  mb,  midway  through  the  depth 
of  the  vortex.  There  is  good  agreement  between  model  and 
observed  winds  at  VBG,  NTD,  and  LIO.  The  forecast  winds  at 
NSI  are  in  error  slightly  as  the  model  has  placed  the  center 
of  the  vortex  near  San  Clemente  Island,  south  of  the  analyzed 
center  near  Santa  Catalina. 

At  0000  UTC  03  May  1990,  the  model  shows  the  sea  breeze 
dominating  flow  in  the  California  bight  at  all  levels  up  to 
875  mbs .  Only  two  reports  are  available  to  confirm  the 
model's  wind  fields,  VBG  and  LIO.  Model  output  winds  are  in 
good  agreement  with  these  observations.   Additionally,  the 
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model  wind  fields  agree  well  with  the  0000  UTC  composite 
Catalina  Eddy  of  Mass  and  Albright  (1989) . 

Figure  64  shows  980  mb  model  wind  fields  and  observations 
at  1200  UTC  03  May  1990.  The  few  observations  available 
cannot  define  the  location  of  the  Eddy's  vortex,  however, 
southerly  flow  at  NSI  at  980  mb  confirms  its  existence 
somewhere  southwest  of  NSI.  The  model  wind  fields  are  very 
reasonable,  and  match  observed  winds  at  VBG,  NTD,  LIO,  and  NSI 
quite  well. 

During  the  period  of  strongest  Catalina  Eddy  activity 
(1200  UTC) ,  the  model  shows  a  region  of  stronger  offshore  flow 
funneled  through  the  passes  between  the  San  Ynez  and  San 
Rafeal  Mountains  at  all  levels  up  to  875  mb.  The  model 
suggests  this  flow  originates  east  of  the  Sierra  Nevada, 
following  the  terrain  as  it  flows  around  higher  barriers  and 
over  the  lower  topographic  features  until  it  reaches  the 
vicinity  of  the  California  bight  (Figure  65) .  The  existence 
of  this  flow  pattern  is  supported  by  observations  at  both  EDW 
and  NID.  Upon  reaching  the  California  basin,  the  wind  maximum 
would  have  considerable  horizontal  shear  (and  associated 
positive  vorticity)  over  the  bight  region.  The  model's  980  mb 
relative  vorticity  field  at  1200  UTC  03  May  (dashed  lines, 
Figure  64)  confirms  that  cyclonic  shear  vorticity  occurs  in 
the  California  bight.  This  vorticity  could  be  a  possible 
maintenance  or  generation  mechanism  for  the  Catalina  Eddy,  and 
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compliment  the  lee  troughing  mechanism  postulated  by  Mass  and 
Albright  (1989) . 
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VI.  CONCLUSIONS 

This  study  determined  that  the  NRL  Limited  Area  Dynamical 
Weather  Prediction  Model  can  be  used  to  study  mesoscale 
atmospheric  phenomena  along  the  west  coast  of  the  United 
States.  In  particular,  the  model  successfully  replicated 
three  instances  of  mesoscale  coastal  phenomena  observed  in 
that  region  during  the  period  0000  UTC  02  May  1990  to  1200  UTC 
03  May  1990.  They  include  land/sea  breezes  as  observed  by 
the  Naval  Postgraduate  School  Doppler  wind  profiler,  a 
Catalina  Eddy  event,  and  its  associated  southerly  surge. 
Although  the  model  had  problems  with  exact  replication  of  near 
surface  low- level  inversions,  in  general  model  wind  and 
temperatures  had  only  small  departures  from  observed  fields. 
Validation  of  the  height  and  moisture  fields  were  inconclu- 
sive, however,  and  further  work  is  required. 
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